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We study theoretically the osmotic pressure of a suspension of charged objects (e.g., colloids, 
polyelectrolytes, clay platelets, etc.) dialyzed against an electrolyte solution using the cell model 
and linear Poisson-Boltzmann (PB) theory. From the volume derivative of the grand potential 
functional of linear theory we obtain two novel expressions for the osmotic pressure in terms of the 
potential- or ion-profiles, neither of which coincides with the expression known from nonlinear PB 
theory, namely, the density of microions at the cell boundary. We show that the range of validity 
of linearization depends strongly on the linearization point and proof that expansion about the 
selfconsistently determined average potential is optimal in several respects. For instance, screening 
inside the suspension is automatically described by the actual ionic strength, resulting in the correct 
asymptotics at high colloid concentration. Together with the analytical solution of the linear PB 
equation for cell models of arbitrary dimension and electrolyte composition explicit and very general 
formulas for the osmotic pressure ensue. A comparison with nonlinear PB theory is provided. Our 
analysis also shows that whether or not linear theory predicts a phase separation depends crucially 
on the precise definition of the pressure, showing that an improper choice could predict an artificial 
phase separation in systems as important as DNA in physiological salt solution. 

PACS numbers: 82.70.Dd, 64.10.+h 



I. INTRODUCTION 

In this paper we study the osmotic pressure of a sus- 
pension of charged colloids or polyelectrolytes in osmotic 
equilibrium with an electrolyte of given composition. Ex- 
amples of such systems abound in our everyday life. They 
occur as dispersion paints, viscosity modifiers, floccu- 
lants, or superabsorbers, to name but a few technologi- 
cal applications jl], ||. They also play a tremendous role 
in molecular biology, since virtually all proteins in ev- 
ery living cell, as well as the DNA molecule itself, are 
charged macromolecules dissolved in salty water |j. A 
great deal of experimental and theoretical research has 
been devoted to their understanding, and several good 
textbooks g], |, |, |, | § and review articles §, [| |<| Q 
are available. 

Arguably the most fundamental thing to know about 
these suspensions is their equation of state, i.e., how 
the (osmotic) pressure depends on other thermodynamic 
variables like macromolecular charge or concentration. 
Within the last hundred years several ingenious ways 
have been conceived for treating this problem on varying 
levels of sophistication. In this article we will be con- 
cerned with Poisson-Boltzmann (PB) theory in combina- 
tion with a cell model approximation for the macroion 
correlations, which we briefly revisit in Sec. II. While 
this does not present the highest level of accuracy or so- 
phistication, it is probably the simplest and up to today 
single most important starting point; it offers a bench- 
mark against which all other theories are compared. In- 
deed, we believe that modern improvements can only be 



fully appreciated once one understands the successes and 
failures of the most fundamental mean field theories. 

Since the nonlinear PB equation can be solved ana- 
lytically only in very few cases p"2[ , its linearized ver- 
sion has always been an important substitute. However, 
the freedom to choose an expansion point and its sub- 
sequent impact on the range of validity and accuracy of 
the linearization has often gone unnoticed. Moreover, 
the computation of thermodynamic properties from the 
ionic profiles computed in linear theory is often based 
on expressions from the nonlinear PB theory or expan- 
sions thereof Jl3], [l^, |l6). This procedure is by no 
means unique and invariably entails internal inconsisten- 
cies. Both these points make it virtually impossible to 
conclude whether any failures of linearized PB theory are 
real deficiencies or avoidable side-effects of a non-optimal 
or inconsistent linearization. 

In this paper we resolve these issues by giving a co- 
herent presentation of linearized PB theory which illu- 
minates the subtle interrelations between the Donnan 
equilibrium, microion screening, linearization and the 
osmotic pressure. In Seolllflwe utilize the functional 
approach to PB theory |17t |18|, 19 and generalize its 
quadratic expansion |2C| ], arriving at a functional which 
yields the PB equation linearized about the electrostatic 
potential value ip. This expansion point will turn out to 
lie at the heart of all those interrelations. For general 
ij) we then derive in Sec. IV an analytical formula for 
the pressure in terms of the ionic profiles. Our novel ex- 
pression replaces the famous boundary density rule |2^| 
from nonlinear PB theory, according to which the os- 
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motic pressure of the suspension is given by the value 
of the microion density at the outer cell boundary. That 
this does not hold in linearized theory may be considered 
as one of the major results of the present work. 

The most common choice of the linearization point ip - 
namely, the potential value in the salt reservoir - suffers 
from several drawbacks, as this value can be very dif- 
ferent from the average electrostatic potential (ip(r)) in 
the colloidal suspension, and it seems more reasonable to 
selfconsistently linearize about the latter |l4, ^2], ^3], p4| . 
We investigate this choice in Sec. [v| and wi 1 prove that 
it is indeed optimal in the sense that (i) charge screen- 
ing rests on the actual ionic strength, such that (ii) the 
crossover between counterion- and salt-screening is natu- 
rally included and (Hi) the limit of large volume fraction 
is correctly reproduced. In brief, all zeroth order effects 
of the Donnan-equilibrium are already incorporated by 
the mere choice of the linearization point, and the lin- 
earized equation now describes higher order effects. We 
will also see that (iv) all other choices of the lineariza- 
tion point overestimate the Donnan effect in lowest order, 
since they violate a rigorous inequality from PB theory 
for the salt content in the colloidal suspension. 

Based on this optimal linearization scheme we then de- 
rive explicit analytical formulas for the osmotic pressure 



of suspensions of charged mesoscopic objects in Sec. VI 



These formulas hold for spherical, cylindrical and planar 
shapes and should thus be useful in a broad variety of 
possible systems. The subsequent Section VII is devoted 



to comparing their predictions with the full nonlinear PB 
theory. 

Within full PB theory the pressure is always positive. 
Whether or not this also holds in the linearized theory 
depends both on the choice of the linearization point ip 
as well as on the precise definition of the pressure it- 
self. We prove that the pressure is always positive for 
symmetric electrolytes if one treats ij) as an independent 
variable. If one does not, the pressure can become neg- 
ative at low volume fractions ^5], ^6], p9[ . 
The implied liquid-gas coexistence - not being present 
on the nonlinear level - is thus clearly an artifact. As a 
striking example we show that even a solution of DNA 
molecules under physiological conditions would be pre- 
dicted to phase separate at all relevant densities. 
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FIG. 1: A solution of charged objects is in osmotic (Donnan-) 
equilibrium with a salt reservoir of given composition. The 
membrane is permeable for small ions only and has to support 
an excess osmotic pressure AP. The mesoscopic objects are 
depicted as spherical (as appropriate for instance for "conven- 
tional" charged colloids or micelles), but our discussion will 
be more general and also apply to cylindrical entities (e.g. 
DNA, actin filaments, TMV viruses) or planar objects (e.g. 
charged membranes, clay particles). 



A. The Donnan equilibrium 

We study a suspension of charged mesoscopic objects 
(henceforth simply referred to as "colloids") dialyzed 
against a salt reservoir of given composition, as illus- 
trated in Fig. [y. This situation is traditionally referred 
to as a "Donnan membrane equilibrium" |30| [nj |32| . 
Much of our discussion will not depend on the shape 
of the colloids, and our final explicit formulas will be 
valid for spherical, cylindrical and planar geometries. 
Even though the microions can traverse the membrane, 
their average concentration differs between the salt reser- 
voir and the colloid compartment, since the latter is al- 
ready occupied by the counterions originating from the 
macroions, which cannot leave the compartment due to 
the constraint of global electroneutrality. This imbalance 
in average densities generates an osmotic pressure differ- 
ence which the membrane has to sustain and which we 
wish to calculate in the following. 

For simplicity we assume the reservoir to be sufficiently 
large such that its ionic strength remains unchanged af- 
ter being brought in contact with the macroion solu- 
tion. This assumption is not necessary, but simplifies 
our discussion of general theoretical issues. How it can 
be avoided is demonstrated in Ref. |3^| . We also note that 
we will not describe the solvent explicitly but rather re- 
place it by a continuum with relative dielectric constant 



II. GENERAL FRAMEWORK 



In this section we start by introducing the physical sit- 
uation we wish to describe - namely, the Donnan equi- 
librium - and its theoretical description in terms of a cell 
model. PB theory is founded on its grand potential func- 
tional, and a brief derivation for the pressure (leading to 
the boundary density rule) is presented. 



B. Cell- model and Poisson-Boltzmann theory 

Theoretical concepts like the cell-model or PB theory 
have long become standard tools, so we will restrict our- 
selves to a brief description and only provide the basic 
equations — essentially in order to introduce our notation 
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and terminology. A recent and more detailed exposition 
can be found in Ref . |33) . 

The cell model approximation attempts to reduce 
the complicated many particle problem of interacting 
charged colloids and microions to an effective one-colloid- 
problem, ft rests on the observation that at not too 
low volume fractions the colloids - due to mutual re- 
pulsion - arrange their positions such that each colloid 
has a region around it which is void from other colloids 
and which looks rather similar for different colloids. In 
other words, the Wigner-Seitz cells around two colloids 
are comparable in shape and volume. One now assumes 
that (i) the total charge within each cell is exactly zero, 
(u) all cells have the same shape and (Hi) for actual 
calculations one may approximate this shape such that 
it matches the symmetry of the colloid (for instance, 
spherical cells around spherical colloids). If the radius 
of the colloids is tq, the cell radius R is chosen such that 
(j) = (r /R) d equals the volume fraction occupied by the 
colloids. Here, d measures the "dimensionality" of the 
colloid in the sense that d = 1, 2 and 3 corresponds to 
planar, cylindrical and spherical colloids. If one (iv) ne- 
glects interactions between different cells, the partition 
function finally factorizes in the macroion coordinates, 
i.e., the thermodynamic potential of the whole suspen- 
sion is equal to the number of cells times the thermody- 
namic potential of one cell. 

Within each cell the small ions assume some inhomo- 
geneous distribution which arises from their interactions 
with themselves as well as with the charged colloid. Com- 
puting the corresponding partition function is imprac- 
ticable, since all ions are strongly correlated with each 
other. Poisson-Boltzmann theory is the mean-held route 
to circumventing precisely this problem. A very powerful 
way to formulate it starts from a thermodynamic poten- 
tial functional belonging to the appropriate ensemble — 
which in our case is the grand canonical one: 



(3VL = 



+ ^n*(r)( log (rn(r)A d ) -1 

i 



(1) 



The meaning of the symbols is as follows: j3 = Xjk^T is 
the inverse thermal energy; ijj( r ) is the local electrostatic 
potential (made dimensionless by multiplication with (3e, 
where e is the positive unit charge); the potential is gen- 
erated by both the fixed charge density ep{(r) (located 
for instance on the colloid surface) as well as the distri- 
butions fii(r) of mobile ions of species i, which have a 
signed valence Vi and a chemical potential pn\ A is the 
thermal de Broglie wavelength of the small ions; and the 
region of integration, V, is understood to be the space 
within one cell that is actually accessible to the small 
ions. The functional minimization of (|lj) subject to the 



constraint of Poisson's equation and charge conservation 
yields the set of Euler-Lagrange equations 



i(r) = A~ d e^' e~ v ^ {r) = 



-Vilp(r) 



(2) 



The rii t b are the concentrations of ions of species i in the 
salt reservoir where the electrostatic potential has been 
assumed to vanish. Combining Eqn. (||) with Poisson's 
equation results in the nonlinear Poisson-Boltzmann dif- 
ferential equation for the potential tp(r) in the region 
within the cell accessible to the microions. After intro- 
ducing the Bjerrum length 1-q :— /3e 2 /e r , it is written 
as 



Aip(r) = -47r l B y^i rii(r) 



(3) 2 Hi V i n i,b c 
— K 



—Vitp(r) 



(3) 



with the Debye screening constant of the reservoir defined 
as 



(4) 



If we reinsert the solution of this equation back into the 
functional (0) and use Eqn. (Q), we obtain its equilib- 
rium (=minimum) value, which is the grand potential of 
nonlinear PB theory: 

f3n cq = ~ [ d d r pi (rty(r) 



d d r ni(r) 



1 + 2 y i ^( r ) 



(5) 



We want to close with the following remarks: The 
above variational principle can be constructed starting 
from the PB equation (see for instance Ref. [^9) and refer- 
ences therein) , but this need not give a unique functional 
p4| and (once f2 has been identified with the grand poten- 
tial) appears like an upside-down explanation for the key 
initial equation (||). However, PB theory can be well jus- 
tified by deriving the functional ([!]) from the underlying 
Hamiltonian. For instance, it can be obtained as the sad- 
dle point of the field theoretic action ]35|, |3^, [37j |38| , as 
a density functional reformulation of the partition func- 
tion combined with a first order cumulant expansion of 
the correlation term pc[| , or from the Gibbs-Bogoljubov 
inequality applied to a trial product state [^3). Those 
approaches also show that PB theory provides an upper 
bound of the exact thermodynamic potential. 



C. The pressure in Poisson-Boltzmann theory 

One advantage of the thermodynamic functional ap- 
proach to PB theory is that it becomes immediately clear 
what the pressure is — in the present ensemble the deriva- 
tive of the grand potential (0) with respect to the volume. 
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It proves convenient to rewrite this in terms of the func- 
tional (Q), which can be achieved as follows. The vari- 
ation upon some change in volume can be decomposed 
into the "orthogonal" changes 



sn 



SV 



Jv 



5m(r), (6) 



where the first term contains any explicit dependence on 
the volume (at fixed ion profiles) and the second part 
is the implicit dependence through the ion profiles (at 
fixed volume). However, since the equilibrium distribu- 
tions make the grand potential functional stationary with 
respect to variations of the density profile at fixed cell ge- 
ometry, the implicit terms vanish. Hence, the pressure is 
just the negative derivative of the grand potential func- 
tional with respect to the cell volume, evaluated at the 
equilibrium profile. The derivative d/dV is understood 
to imply a movement of the outer neutral cell boundary, 
which shall be located at r = R, and which only oc- 
curs in the boundaries of the volume integral in Eqn. (Q). 
Note that rewriting the electrostatic energy in terms of 
the densities gives a double integral, and the product 
rule then cancels the prefactor 1/2 in front of the term 
describing the electrostatic energy. Putting everything 
together, one arrives at 



(3P 



dV 



dV 



(7) 



i 

= 5>, b e"^) 

i 



v^(R)+\og (rn{R)A d ) -1-in 



cq 

(8) 
(9) 



This is the well known result [plj that within PB theory 
the pressure is given by the sum of the ionic densities at 
the cell boundary. It actually holds beyond the mean- 
field approximation |]39f , but this will not be our concern 
in the following. Its generalization for more complicated 
cells can for instance be found in Refs. jll| p3| . 

Eqns. (||,^|) give the pressure acting within the 
macroion compartment. The excess osmotic pressure 
across the membrane is the difference between this pres- 
sure and the pressure in the salt reservoir. The latter 
is obtained with comparative ease, since the electrolyte 
is homogeneous and hence the route via a density func- 
tional is unnecessary. On the same level of approximation 
as above it is given by the van't Hoff equation 



nib- 



(10) 



This implies for the excess osmotic pressure across the 
membrane 



f3AP = ^ni, 



b(e 



-Viip(R) 



-1) > 0. (11) 



The last inequality follows from e x > 1 + x and the fact 
that the salt reservoir is neutral, i.e., Yli v i n i,b = 0. 
Hence, within PB theory the excess osmotic pressure is 
always nonnegative. 



III. THE GRAND POTENTIAL IN 
LINEARIZED THEORY 

When studying linearized PB theory we want to ben- 
efit from the same thermodynamic coherence as in the 
nonlinear case. We can achieve this aim by likewise 
founding the key equations on a suitable grand poten- 
tial functional — a possibility that has previously been 
pointed out by Lowen etal. p(| . 

First observe that a functional leading to a linearized 
version of the PB equation necessarily is quadratic in the 
densities. The only part in Eqn. (Q) for which this is not 
already true are the entropy terms 7Jj(r)[log(ni(r)A d ) — 
1]. Let us now introduce a set of densities hi by 
"Boltzmann-weighting" the reservoir densities n^b with 
an as yet unspecified potential value ip: 



-Vilp 



(12) 



If we expand the entropic terms in Eqn. (|lj) about the 
points hi up to quadratic order we obtain 



1 + 2(1 + ^) 



rxi(r) /n*(r) 



n, 



V hi ) 



(13) 



We will refer to this expression as the grand potential 
functional of linearized PB theory, since its functional 
minimization (again, under the constraint of Poisson's 
equation and charge neutrality) leads to 



n t (r) = n^e-** [l - v^r) - . 



(14) 



This is the Boltzmann relation (g) linearized about the 
potential value ijj. We want to stress right from the be- 
ginning that we leave the value of tp unspecified for the 
time being. Linearized PB theory is not unique, it is 
a one-parameter family labelled by the expansion point. 
The by far most common choice found in the literature is 
■ijj = (strictly speaking, the potential in the salt reser- 
voir), but this is not the only conceivable (let alone op- 
timal) possibility. In Sec. we will come back to this 
issue in greater detail. 

Combination of Eqn. ( |l4]) with Poisson's equation 
yields the linearized Poisson-Boltzmann equation 



Atp(r) 



(14) 



. E l v i n i,b 

k 2 (tp(r) - ipi), 



E 4 v 1 n i,b 



-(ip{r) - ip) 



(15) 
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with the renormalized screening constant R and the in- 
homogeneous term ip{ of the differential equation defined 

as 



and ipi = ijj + 



,2 T,i v i n i 

E,i v i n i,h 

J2i V i fn 



4art B vf «i (16) 
(17) 



Note in particular that k appears as a screening constant 
calculated with the Boltzmann- weighted densities, it 
is hence different from the secreening constant in the salt 
reservoir. For the special case of a v : v electrolyte this 
simplifies to 



Aip(r) "= K 2 ip(r) — n z ib tanh(w-0) 

v 



with 



-2 v:v 2 



k cosh(t>?/;) > k 



(18) 



(19) 



If we reinsert the solution of Eqn. (|T^) into the functional 
( |l3| ) and use Eqn. ( |l4| ) we obtain the equilibrium grand 
potential of linearized PB theory: 



/9fi 



lin.cq 



V 



+ 7 5^«ini(l + «i$)(V(r)-$). (20) 



The angular brackets 



denote the spatial average 



to the ions. 



over the part of the cell volume accessible 



IV. THE PRESSURE IN LINEARIZED THEORY 

As mentioned in the introduction, the computation of 
the pressure in linearized PB theory has often been based 
on formulas origina ting from nonlinear case or expansions 
thereof |l| [lj, |l5| , [Tq |. For instance, one could use the 
predictions for boundary potential or density from lin- 
earized theory and insert them into formulas (JsJ) or (|j), 
respectively. However, although both formulas coincide 
on the nonlinear level, they yield different results once 
the linearized equation is used to compute ip(R) or n,(P), 
since Eqn. (|J) no longer holds. 

Here we circumvent this source of inconsistency by 
avoiding any recourse to resu lts from nonlinear PB the- 
ory. Analogously to Sec. II C we base the pressure on the 



volume-derivative of the grand potential of linearized PB 
theory, Eqn. (^0|). This leads to a formula which gives 
the pressure as a function of the ionic profiles and which 
replaces the boundary density rule (0). 



A. Relevant thermodynamic variables 

Before we differentiate the grand potential we would 
like to pause for a moment and discuss the issue of 
relevant thermodynamic variables, since the pressure- 
formula will turn out to depend on our choice for them. 

The grand potential of nonlinear PB theory depends 
on volume V, temperature T, colloid charge Q, volume 
fraction and the set of chemical potentials In ad- 
dition to these variables the grand potential underlying 
linear theory also depends on the linearization point, ip. 
This proves a relevant issue since ip itself may depend on 
volume, see the following Sec. [v|. For reasons which will 
become clear later it turns out to be extremely useful to 
consider ib as an independent variable rather than insert 
the functional dependence if>(V) into the potential and 
thereby eliminate ib. 

With this in mind, the thermodynamic definition of 
the pressure within linearized PB theory becomes 



P 



(i) 



lin 



/ anim,cq(v,r, 
^ dV 



A) 



(21) 



i. e., the partial derivative of the potential (|20| ) with re- 
spect to the volume, keeping all other variables fixed. 
This amounts to the following procedure: If one wishes to 
calculate the pressure at a given volume, one first chooses 
the desired linearization point ip(V) at this volume, but 
fixes it subsequently. Then one measures the change in 
the grand potential upon slightly changing the volume. 

However, one could also argue that ^im,eq(^ • ■ ■ ) : = 
^iin,eq(^ • ■ ■ j ^00) is the desired grand potential. In 
this case ip is not regarded as an independent variable, 
but is removed from the description by substitution. The 
pressure from the derivative of this potential is then 



P, 



(2) 



aniin,cq(y,r,---) \ 
ov j 

V / T, ... 

'dniin.cqCv,? 1 ,... ,^{V)) 

dV 



lin 



V.T,. 



(22) 

drp(V) 
dV ■ 



This is the total derivative of the grand potential 
^iin,eq(y • ■ ■ , ip{V)) with respect to volume, and it dif- 
fers from by an additional term which stems from 
the volume dependence of ip. The incorporation of the 
constraint ip(V) amounts to the restriction of possible 
(V, tf>) values to a one-dimensional submanifold. One 
hence compares values of the grand potential at neigh- 
boring points on that submanifold, along which one dif- 
ferentiates. Fig. H gives a schematic illustration of the 
difference between these two definitions. 

After these considerations we can now proceed with 
our aim to derive a formula which gives the pressure in 
terms of the solution of the linear equation. 
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SXVMV)) = £i(V) 




FIG. 2: Relation between the two pressure definitions from 
Eqns. @ and ( p2j ) discussed in the text. In the first scheme 
the derivative is taken in the direction of increasing V along 
a coordinate line of constant tp; in the second scheme the 
derivative is taken tangential to the submanifold Q(V, tp(V)), 
which emerges from restricting the values (V, tp) to the ones 
given by the defining equation tp(V) of the linearization point. 



and one readily obtains using Eqns. (H) and (|l4|) 



on 



lin,cq 



lin 



dip 



dip 



cq 



vfrn. 



(25) 



mula for Pi? : 



Combining this with Eqn. (123) results in an explicit for- 



(i) 



2 dlA v ^ ' 



v? m 



(26) 



Note that unlike Pj£ this expression depends on the 
whole potential distribution, and not just on the bound- 
ary potential or the boundary densities. 
Equations 



paper. They rep 



1Q and (26) are a first key result of this 
ace the boundary density rule (^|) from 
nonlinear PB theory whose validity in the linear case 
has often falsely been taken for granted. In the second 
part of this article we will use these equations to derive 
explicit formulas — once we have established an optimal 
linearization point tp, which will be the topic of the next 
section. 



B. The derivative of the functional 

The same line of reasoning that led to Eqn. (R) can 
be employed to rewrite the derivative of the equilibrium 
grand potential as a derivative of the grand potential 
functional. If one again remembers that d/dV just corre- 
sponds to a movement of the outer neutral cell boundary, 
the pressure definition (|2l|) and Eqn. (|l4|) give 



[-iH ni{R)ip(R) 



y 



+ (l + v^n^R) 
$> 4 , b e-^ 

i 

i 



1 (n z (R)Y 



2 fii 
l-Vi(ip(R) 



cq 



-vf(i,(R)-i,y 



i 



7ii(R) rij 



ni(R) 



(23) 
(24) 



Eqn. (23) replaces Eqn. (|3|) and is easily recognized as 
its quadratic expansion about tp. Similarly, Eqn. (24) 
replaces Eqn. ([)[). Since x + l/x > 2, the expression (24) 
is larger than ^2ini(R), i.e., the expression that would 
follow if the boundary density rule (|^) would also hold in 
the linear case. 



In the case of pressure definition (22) we need the ex- 
plicit volume dependence of Oii n oq , which can again be 
rewritten as the derivative of functional r2n n , Eqn. (13). 
From Eqn. ( p^ ) follows immediately dTii/dip = —Vifii, 



V. OPTIMAL CHOICE OF V 

Up to now we have not specified the linearization point 
tp; rather, we have emphasized that its choice is largely at 
ones disposition. However, not all choices may be equally 
successful. In fact, the range of validity of linearization 
depends strongly on the choice of ip, since it is the dif- 
ference between the potential and its linearization point 
that is required to be small and not the potential itself. 
This is particularly important in concentrated suspen- 
sions where this difference can be small even when the 
total potential is quite large. 

In this section we identify an optimal linearization 
scheme — which however first requires a clarification of 
what "optimal" is supposed to mean. It proves instruc- 
tive to first study the obvious (but futile) attempt to 
base optimality on a min imiz ation of the grand poten- 
tial, as we discuss in Sec. V A . This sheds some surpris- 
ing light onto traditional linearization. 



VB 



we 



In Sec 

will then discuss a scheme based on the selfconsistcntly 
determined a vera ge potential and present its optimal as- 
pects in Sec. V_C, which even though this approach has 
been used in the past |l4|, plL ^3[ |24| have largely gone 
unnoticed. 



A. Minimizing the grand potential? 

One may try to obtain an optimal expansion point by 
looking for the value of ip which minimizes the grand po- 
tential of linear theory. Let us thus set the derivative 
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of fiiin.eq, Eqn. (f20|) , with respect to ip to zero. Us- 



ing Eqn. ( |2q ) and remembering that {{ip(r) — ^) 2 ) is 
nonzero unless the profile is completely flat, we see that 
the value of ij) at the extremum is given by the solu- 
tion of vfrii t b e~ Vi ^ = 0. The left hand side is a sum 
of strictly monotonically decreasing functions. If ions 
of both sign are present, the left hand side hence falls 
monotonically from +00 to —00 and the equation has a 
unique solution. Unfortunately, however, this monotonic 
decrease of the derivative also implies that this solution 
corresponds to a maximum of the grand potential. For 
a symmetric v : v electrolyte the solution is for instance 
given by ip = 0. The most widely used choice of the ex- 
pansion point hence gives the largest grand potential of 
all possible linearization schemes (this can for instance 
be verified in Fig. 1 of Ref. p^j , which shows the ther- 
modynamic potential for various linearization schemes). 

We hasten to remark that the above finding has to 
be put in the correct perspective. Minimization is only a 
meaningful venture if one can be sure about the existence 
of a lower bound |57|]. For instance, the PB functional 
is bounded below by the exact thermodynamic potential 
of the restricted primitive model. Minimizing the func- 
tional stems from the desire to get as close to this result 
as is possible within a mean field description Q. For 
the parameter ip from linearized PB theory such a lower 
bound for the functional cannot be constructed, hence 
minimization is meaningless. And even if there were a 
bound, there would be no reason to approach it — unless 
one knows that it is favorably related to the actual ther- 
modynamic potential. 



B. Expansion about the average potential 

Having seen that a minimization condition on fin n ,eq 
is not successful, we will approach the problem from a 
different direction. If we average Eqn. (UJ) over the cell 
volume and use Eqn. (112), we find 



(nj(r)) = m - Vim(ip(r) -tp). 



(27) 



If we were to choose (ip(r) — ■ip) = 0, the second term 
would vanish and the expansion points for the densities, 
fij, would coincide with the averages (rii(r)}. It clearly 
makes sense to expand about these average values, since 
then the differences between actual value and expansion 
point can be kept small throughout the cell, i.e., in the 
whole region in which linearization must work. Hence, 
the choice 



t/j = t/j, 



opt 



(28) 



is a particularly suitable one, which in anticipation of 
the results from Sec. V C we have labelled with the index 
"opt" . 

At first sight this particular choice may seem difficult 
to work with, since the average (VK r )) determines the 
value of V^opt, but the potential ip{r) to be averaged is in 



turn the solution of the equation linearized about V'opt- 
However, the necessary selfconsistency is readily fulfilled, 
since (i^(r)) (and thus V> op t) does only depend on the 
state point and not on the actual ionic profile: For rea- 
sons of electroneutrality the total charge of all ions in the 
cell must be the negative of the colloid charge, and so we 
have 



Q 
V 



Vi{n,i 



(r)> 



i n»,b 1 



-l>zl/>o 



(29) 



Satisfying this equation is necessary and sufficient for the 
validity of Eqn. (p8|), so we can determine V> pt also from 
Eqn. (E9). If the electrolyte contains ions of both sign, 
the right hand side monotonically decreases from +00 to 
—00 as a function of V'opt, so a solution always exists and 
is unique. For the special case of a v : v electrolyte this 
equation reads 



Q/vV 
2n h 



sinh^Vopt)- 



(30) 



C. Optimality of ifj op 



The linearization scheme from the above Sec. VB has 
been put forward several times in the past O, |22|, |23|. p4| . 
but various of its special properties have gone unnoticed. 
In this section we show, in which sense this scheme can 
be regarded as optimal. 

Since the choice of tp op t implies that the expansion 
points fii coincide with the average microion densities 
(nj(r)), this translates to the corresponding screening 
constant K opt which - as Eqn. ( |l6| ) shows - is now also 
calculated with the average ionic strength within the cell: 



^opt = 4^B^>- 



i(r)>. 



(31) 



Due to the Donnan effect the latter is different from the 
one in the reservoir, and it is more appropriate to have 
a linearization scheme which derives its screening con- 
stant from the actual ion densities in the macroion com- 
partment. In fact, describing the colloidal system in an 
integral-equation approach together with the MSA clo- 
sure yields an effective pair potential between colloids 
that is a screened Coulomb potential with the screening 
constant being equal to /c pt — see for instance Ref. O] . 

Let us make the above remarks more explicit by spe- 
cializing to a symmetric v : v salt of concentration 
rib- The average counterion concentration is denoted by 
n c := \Q\/vV. In this case the optimal screening constant 
can be expressed as 

-2 (1?) .2 



„ opt - K COsh^opt 

(§) 



8ne B v 2 n h ^l + {n c /2n h ) 2 
4ttI b u 2 V^1 + ( 2 «b) 2 . 



(32) 
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This expression differs from the screening constant of the 
salt reservoir by also incorporating the contribution due 
to the counterions. It is illuminating to study its limit- 
ing behaviors at high and low salt. In the salt dominated 
regime, i. e. n c <C 2nb, the screening constant K opt is 
practically identical to the reservoir screening constant n. 
This is the limit in which the Donnan effect is strongly 
diminished, and Eqn. ( |30| ) also shows that in this case 
"0opt — * 0. In the opposite limit of counterion domina- 
tion, n c 3> 2rib, the screening constant K opt is essentially 
determined by the average counterion concentration, and 
it will be substantially larger than the screening constant 
of the reservoir. In this limit the Donnan effect is strong 
and i/'opt is significantly different from the zero. The 
traditional linearization scheme ip = will fail to ade- 
quately describe this case, even though the ionic profiles 
can be very flat and thus amenable to linearization. In 
summary, we see that linearizing the PB equation about 
'tjj = invariably ignores the counterion contribution to 
the screening, whereas the expansion about the selfcon- 
sistently determined average potential includes the coun- 
terions in a way that yields the correct behavior in the 
limits of high and low salt. 

The fact that a strong Donnan effect goes along with 
a value of tp op t substantially different from the reser- 
voir potential is also to be expected on more fundamen- 
tal grounds: The difference in microion concentrations 
across the membrane characteristic of the Donnan equi- 
librium implies a concomitant jump of the electrostatic 
potential, referred to as the "Donnan-potential" . Its 
magnitude can likewise be used to measure the strength 
of the Donnan effect. But a large Donnan-potential also 
results in a large average potential. This close connec- 
tion between the Donnan-equilibrium and the optimal 
linearization scheme lead the authors of Ref. j2J] to sug- 
gest the name "Donnan-linearization" |5£| for the scheme 
from Eqn. @. 

The above findings can be succinctly reformulated in 
the following way: By choosing the optimal linearization 
point ipopt the Donnan-equilibrium is automatically de- 
scribed correctly to lowest order. The solution of the 
linearized PB equation then provides the higher order 
corrections. We now show that every other lineariza- 
tion point gives a description that violates an important 
zcroth-order inequality valid on the nonlinear PB level. 
Let us first specialize to a v : v electrolyte. If we denote 
by n c the average density of counterions and by n s the 
average density of salt molecules within the cell, we have 
within Poisson-Boltzmann theory 



find 



(n c 



n s )n s ( = } nl[l-v(ip(r)-tp}] [l + v(ip(r) 
= n 2 h [l-v 2 (Mr)~4>) 2 ] < nl 
>o 



(34) 



This is the inequality ( |33| ) but with reversed inequality 
signl Observe that linearizing about ijj op t is the only 



scheme that does not violate (p3[), since Eqn. ( |28| ) ren- 
ders ([34]) an equality. The importance of these two in- 
equalities can be further unveiled by solving them for 
the average salt concentration n s inside the suspension, 
which gives the sequence 



IPB 



> n s 



LPB, Vo 



> n s 



LPB 



with 



2n s 



LPB, -0o 



(35) 



(36) 



In other words: Linearized PB theory gives a lower bound 
to the average salt density within the cell, and ^opt gives 
the largest lower bound. Since underestimating n s means 
overestimating the Donnan effect, the ^ pt-scheme gives 
as close a representation of the Donnan equilibrium as 
is possible in a linear theory. Observe also that the 
prediction ( |36| ) for n s from Donnan linearization is in- 
deed the well known formula describing the Donnan ef- 
fect after neglecting activity coefficients (see for instance 
Ref. |tt],|4jj). In that sense the set of inequalities ( |35| ) lies 
at the heart of the Donnan equilibrium and distinguishes 
tpo-pt as optimal. 

It is possible - although not straightforward - to ex- 
tend the above considerations to general electrolyte com- 
positions. In that case one has to look at more refined 
combinations of densities, also taking into account to 
what percentage a particular species is represented in the 
ionic mixture. One possibility is the following: Define the 
average 

> 1 ( n i( r )) nib /Q7 ^ 

ai := 2^ftlog— with pi := — , (37) 



which is inspired by a similar procedure employed when 
defining mean activity coefficients |h|. Within PB the- 
ory we can readily derive 



(n c + n s )n s = (n+(r)) (n_ (r)) 



(2) 



-vi/i(r)\ /„«i/>(r)\ 



> nt e 



' b ■ 



(33) 



where Jensen's inequality |39[ has been used. Eqn. ( |33| ) 
states that the mean ionic molarity |Io| within the sus- 
pension is greater than in the reservoir. On the other 
hand, within linearized Poisson-Boltzmann theory we 



i 

> -£>«<Wr)> = 0, (38) 



where Jensen's inequality has again been used. The last 
step follows from the charge neutrality of the salt reser- 
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voir. However, within linearized PB theory we get 



i 
i 

< -^JHVityW-j)) = 0. (39) 



Here, the elementary inequality log(l + x) < x has been 
used. It is straightforward to see that the relations d38| ) 
and ([39]) reduce to ( |33"1 ) and (|34| ) for a v : v electrolyte. 
Nonlinear and linearized PB theory again lead to con- 
flicting inequalities in all but one case: If ip op t is chosen 
as the linearization point, the PB inequality is not in- 
fringed. 



VI. EXPLICIT EXPRESSIONS AND 
APPROXIMATIONS FOR d-DIMENSIONAL 
CELL MODELS 



In Sec. IV we gave expressions of the pressure in terms 
of the charge- or potential-profiles that were based on the 
derivative of the grand potential. In this section we will 
substantiate these results by inserting the actual solu- 
tion of the linearized PB equation for a d-dimensional cell 
model |Q . We will place particular emphasis on the for- 
mulas that emerge from using the optimal linearization 
point V'opt that we have discussed in the last section. 



A. Analytical formulas for the pressure 

In the Appendix we outline how the linearized PB 
equation can be solved for a general c£-dimensional cell 
model (where d — 1, 2, and 3 corresponds to planar, 
cylindrical and spherical macroions, respectively). If we 
insert the final expression ( A5) into the pressure equation 
(E||), we get the explicit formula 



lin 



V > .1 Vi m / 



where the variable 7 is defined as 

l-4> 



(40) 



(41) 
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and D is given in Eqn. ( 

linearization —Q/V = Ei v i^i according to Eqn. (29) 
such that Eqn. (Effl) further reduces to 



dDy/4> 

In the case of Donnan 



(42) 



The case of the alternative pressure definition (|22|) is 
less straightforward, since this depends explicitly on the 
whole potential distribution tp(r) and not just just on its 
boundary value. According to Eqn. (|2^) the hard part is 
the integral over (ip(v)— V> pt) 2 , which is very unwieldy, as 
it contains products of modified Bessel functions. How- 
ever, one can always numerically integrate this expression 
or - alternatively - numerically differentiate the grand 
potential fiii n ,eq from Eqn. (|20|). 

The explicit pressure formulas (ff"o| ) and ( f42|) are a fur- 
ther key result of this paper, and in the remainder we will 
study some of their consequences. Unfortunately, due to 
the algebraic complexity their properties cannot be read- 
ily seen. In the next subsection we will therefore spend 
some time to study their analytical behavior in a few im- 
portant limiting cases. Finally we will provide graphical 
illustrations for the pressure formulas (fh)|) and (^) as 
well as the p ressu re definition ( p2| ) for several exemplary 
cases in Sec. VII. 



B. 



Pressure bounds, limiting behavior and 
expansions 



The pressure formula ( p3|) is quadratic in ip(R), and it 
is easily checked that it takes its minimum value at the 
inhomogeneity if)\ defined in Eqn. ([H]). From Eqn. (A5) 
and the definition of 7 in Eqn. ( fTi] ) we see that this for- 
mally corresponds to 7 = 0, and indeed the pressure 
expression (40) attains its minimum value there. For a 
symmetric v : v electrolyte this implies for the excess 
osmotic pressure 



2n h 



B(P {1) 

^ \ lin, min 



2n h 

— cosh(vip) — — sinh(u^) tanh(u^) — 1 



cosh(w^) — l] 
2 cosh(w-0) 



> 0. 



(43) 



Hence, in the symmetric case the pressure from lin- 
earized PB theory defined via Eqn. ( pi"] ) is nonnegative 
for every chosen linearization point ip. However, this 
does not hold for general electrolytes. It may be veri- 
fied that for the asymmetric case of a 1 : 2 electrolyte 
the excess pressure is negative if if) has the same sign 
as the divalent ion species and |$| € [0; 0.6309]. In 
Donnan-linearization this comes down to the requirement 
that Q has the same sign as the divalent species and 
< \Q\ < 1.596nt,V^, where rib is the reservoir density 
of the monovalent species. Essentially, the pressure can 
become negative when the counterion content within the 
cell is overwhelmed by the salt ions. 

Let us now restrict again to a symmetric v : v elec- 
trolyte and define the parameter 



\Q\/vV 

2n h 



2n h ' 



(44) 
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This ratio between average counterion concentration and 
salt reservoir concentration measures whether we are in 
the salt dominated regime (9 small) or the counterion 
dominated regime (9 large). Using 9, Eqn. (|32|) can be 
rewritten as K„ pt = k 2 V1 + 9 2 . The definition ( p4[ ) shows 
that 9 is small if either V is large (i, e., the volume frac- 
tion of colloids is low) or n D is large (i.e., much salt has 
been added to the system). In both cases R opt R becomes 
large and we may exploit the asymptotic behavior of the 
modified Bessel functions jl2j to approximate the quan- 



tity D from Eqn. ( A4 ) according to 



D 



— ^d/2(n ptr Q )l d / 2 (K opt R) 
„ K d / 2 (^ Qpt r ) pRoptfl 

y/2-KKoptR 

tr >l e *o*(R-ro) 



Hence, 7 oc 1/D vanishes exponentially. The pressure 
P,^* will thus approach its minimum value computed 
above. Since in this limit ^ Q pt — * 0, we may expand 



Eqn. (43) for small tp pt and obtain 



f3AR 



(i) 

lin.min 



(^opt) 4 <g> ^ 



2n h 



(45) 



showing that the excess osmotic pressure (measured in 
units of the reservoir concentration) vanishes as the 
fourth power of 9. We note in passing that the lowest or- 
der calculation based merely on the average counterion 
concentration n c from Eqn. ( |36] ) (or, alternatively, the 
boundary density from salt-free PB theory jl6| ) gives in- 
stead the asymptotic behavior 9 2 /2, i. e., only a quadratic 
dependence. 

The full expression ( |42"| ) contains a term which appears 
to be of second order in 9, namely 9 2 j 2 /2. Since 7 van- 
ishes exponentially, one could be led to the incorrect con- 
clusion (by "prematurely" terminating the expansion at 
this point) that the pressure in linearized PB theory also 
vanishes exponentially jlj| . While it is in fact true that in 
the full nonlinear theory the pressure vanishes exponen- 
tially in the limit of salt-domination (a simple plausibil- 
ity argument for this can be found in Ref. p3[), this does 
not hold for the linearized theory [plj , where the expo- 
nential decay of the term 9 2r y 2 /2 is masked by the quartic 
asymptotic from Eqn. (]45|). In any case, one has to be a 
little bit careful about the physical meaning of this limit. 
It indeed correctly describes a single charged colloid im- 
mersed in an electrolyte. However, a dilute suspension of 
many colloids is not well described in this limit, since the 
lack of mutual repulsions renders some basic assumptions 
about the cells questionable. Furthermore, one must keep 
in mind that our formulas only give the osmotic pressure 
of the suspension due to the microions. Even though this 
contribution most often dominates just because there are 
many more microions than macroions, the contribution 
of the latter has to become dominant once the microion 
term vanishes exponentially. 



The limit of large volume fraction, — > 1, requires a 
little more care. Using R — rQ(p~ 1 l d , an expansion of D 
around 0=1 gives to lowest order 



D — y K^+i^o) I„(x ) + I„+2(a;o) 



+ l u+1 (x ) [K„(x ) + K„ +2 (x )J I x (cT 1 ^ - 1), 

with v = d/2 — 1 and xq — Wopt^o = «To(l + 2 ) 1 / 4 . 
Since 9 oc 1/(1 — 0), both 9 and xo diverge as (j> — > 1. 
Xq/2 times the expression in curly brackets approaches 
1 - c/xl with c = (d 2 - l)(d 2 + 3)/64 §2). The term 7 
therefore behaves asymptotically like 



1 



7 



d(l 



c/x 2 )(^ - 1)V$ 



1 + 



(46) 



This shows that 1 — 7 2 is of order 1/9. Using this, it 
is now straightforward to show that to lowest order the 
high volume fraction limit of Eqn. ( f4(i|) is given by 



0P& 



(47) 



This equation is simple to interpret: It states that 
the pressure is given by the average counterion density. 
This is reasonable, because in this limit the ionic pro- 
files become flat and thus the electrolyte ideal. Indeed, 
Eqn. (||^) merely states that the osmotic coefficient goes 
to 1. We hasten to add that Eqn. (^) only demonstrates 
the proper behavior within the cell model. Real suspen- 
sions crystallize at large enough volume fraction, a tran- 
sition which can only be described correctly once one 
accounts for the ordered phase as well — see for instance 
Refs. §§§§§§. But in order to get the 
phase boundary right, one of course also needs a good 
estimate of the grand potential in the fluid phase, so the 
correct scaling of the optimal linearization at high vol- 
ume fraction is after all practically important. 

Within Donnan-linearization positivity of the pressure 
does not generally hold for the second pressure definition 
(p2|), not even for symmetric electrolytes. For a v : v 
electrolyte this can be seen by differentiating Eqn. ( |30| ) 
with respect to V and inserting into Eqn. (|26]), whereby 
one obtains 



^VMin Mill 

2n h 



2 sinh 2 vijj, 



opt 



2coshi/; pt 



<(^(r) 



■0opt) ), 

(48) 



which is clearly negative. In the limit of low volume frac- 
tion or high salt ^ op t tends to zero, and expression fl4S|) 
scales like (t>V> pt) 2 = 9 2 times factors that do not tend 

to zero. It hence does not vanish as quickly as pfr" ~ 9 4 
(2) 

and the pressure P,- must become negative in this limit. 

(2) 

We finally note that for asymmetric electrolytes P}. ' can 
also become negative, but will not always be smaller than 



P, 



(i) 



lin 
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Let us close with a few remarks on the pressure in the 
traditional linearization scheme ip = 0. Since the lin- 
earization point is volume-independent, the two pressure 
definitions (^l|) and (|2^) coincide. If one calculates the 
excess osmotic pressure by inserting -0 = into Eqn. (23), 
one finds that the zeroth order is canceled by the reservoir 
pressure, while the first order drops out due to reservoir 
electroneutrality, giving 



(A5) (Q/Vf 



2 E* vfn ith 



7 



(49) 



In the special case of a v : v electrolyte this simplifies 
further to /3AP\- m /2nb = 9 2r y 2 /2 (recognize the "second 
order" term from above!). Eqn. ( (49|) shows that this pres- 
sure is always positive. Moreover, in the limit of low <p or 
large salt it vanishes exponentially as 7 2 ~ exp{— 2kR} 
- just as in full PB theory. Even though O pt also ap- 
proaches zero in this limit, the corresponding pressure 
asymptotics is dominated by the way in which it does so, 
yielding a power law as discussed above. In the opposite 
limit of large <p we have seen above that 7 — > 1, hence 
for the v : v case the pressure using ip = O-lincarization 
asymptotically behaves like (3APn n ~ n^/4rib, i. e., it di- 
verges quadratically in n c and not linearly as it should. 
We finally note that a quadratic expression like (49) has 
recently been rederived in Ref. 47 by using a pertur- 
bative expansion of the PB equation and restricting to 
lowest order. However, in this approach it is difficult to 
see that this is in fact thermodynamically consistent. 



VII. EXEMPLARY COMPARISON WITH 
NONLINEAR PB THEORY 

In this section we provide a comparison of the above 
pressure formulas with the results from nonlinear PB the- 
ory |52], which is intended to clarify and illustrate the 
findings from the last sections. Apart from the pres- 
sure, we will also calculate the compressibility, which is 
defined as = y (dV/ dAP) . For a graphical represen- 
tation it is however more convenient to plot the reduced 
inverse compressibility d(3AP/dn c — l/n c k^TJ(f , where 
the average density of counterions is again defined by 
fi c = \Q\/vV. 

In order to underline the generality of the formulas 
derived above, we will use the following three different 
situations: (i) spherical colloids in a 1:1 electrolyte, (n) 
cylindrical colloids in a 1:1 electrolyte, and (Hi) spherical 
colloids in a 1:2 electrolyte. 

The first system we study consists of spherical colloids 
of charge Q — 3500 e and radius ro = 133 nm immersed 
in an aqueous solution (i.e., £b = 0.714 nm) which is 
dialyzed against a 1:1 electrolyte of rather low molar- 
ity 3.6 /jM. Motivated by the prediction of a gas-liquid 



phase separation at these parameters , the authors of 
Ref. p3| used them as an illustration which pushes lin- 
earized theory to the limits of validity, making the way 
in which it deviates from PB theory particularly visible. 
Lowering the colloidal charge entails a successively bet- 
ter agreement between linear and nonlinear theory, as 
can also be seen in Ref. 

Fig. H shows the pressure (left) and the inverse reduced 
compressibility (right) as a function of volume fraction 
(p. The solid line corresponds to the solution of the non- 
linear PB equation. The pressure is given by Eqn. (H) 
and the excess pressure is always positive (see Eqn.(|ll|)). 
The same holds for the compressibility J63[ . The dashed 
and dotted curves correspond to the solutions from lin- 
earized PB theory. For -0 = both pressure definitions 
( pl| ) and ( p2| ) coincide, but in this case the predictions 
are clearly off from the PB result except at very low vol- 
ume fraction. In contrast, using the optimal linearization 
point op t from Sec. |VB| brings about the correct behav- 
ior at high volume fraction, i.e., in the regime where 
the cell model is particularly appropriate and lineariza- 
tion should indeed work because of relatively flat ionic 
profiles. 

The pressure definition (|2^) - the total derivative of 
the grand potential fi(V, T, . . . , 4> op t{V)) with respect to 
volume - corresponds to the one employed in Ref. |24| . It 
can be seen to lead to negative pressures and compress- 
ibilities at moderately low volume fractions, which would 
imply a segregation into a dense and a dilute colloidal 
phase. Such a phase transition has also been claimed in 
other recent theoretical works [^5[ ^7[ ^8) which likewise 
are essentially based on linearized PB theory. However, 
the similarity of the resulting phase diagrams, the ap- 
pearance of essentially the same "volume terms" in the 
grand potential (which here are responsible for the ef- 
fect), and the fact that full PB theory does not show 
any phase transition, lead the authors of Ref. 24 to 
question these claims of a gas-liquid coexistence in such 
systems as a spurious side effect of linearization. In a 
similar spirit, Diehl etal. |29|| show (within a general- 
ized Debye-Huckel-Bjerrum approach (4^, [l9|) that ex- 
plicitly re-incorporating the effects of counterion conden- 
sation, which are neglected when doing the linearization, 
removes (or at least strongly suppresses) the phase tran- 
sition otherwise clearly visible. We remark that Fig. 1 
of Ref. ^9|, showing the excess osmotic pressure as a 
function of volume fraction for a salt-free suspension of 
spherical colloids with varying bare charge, can be re- 
produced almost quantitatively within the cell model by 
using Donnan-linearization and employing the pressure 
definition (|22]). While these warnings about the dangers 
of linearization are certainly well made, we want to point 
out that things are in fact even a little more subtle: Not 
even linearized theory needs to show the phase transi- 
tion. Whether or not it does so depends crucially on the 
pressure definition as well as the linearization point. For 
ip = there is no phase transition. For tp — V-'opt there 
is a phase transition only if pressure definition (E2l) is 
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FIG. 3: Pressure (left) and inverse reduced compressibility (right) as a function of volume fraction cj> for a solution of spherical 
colloids having a bare charge Q = 3500 e and a radius ro = 133 nm. The system is in contact with a dilute 1:1 electrolyte of 
concentration 3.6 /iM and the Bjerrum length is £b ~ 0.714 nm, corresponding to water at room temperature. The solid curve 
is the result from nonlinear PB theory, the bold dashed curve is the first pressure definition ( p3| ) combined with the optimal 
linearization from Sec. |v] (leading to the explicit expression (|42|)) while the dotted curve combines this scheme with the second 
pressure definition (^). The fine dashed curve uses ip = 0, in which case both pressure definitions coincide, leading to the 
explicit formula (EsJ). 



used. The definition (|2l| ) - which is based on the partial 
derivative of the grand potential with respect to volume 
- yields a positive compressibility. 

Let us repeat that this pressure definition ( pl| ) could 
only be written down once the linearization point has 
been recognized as an independent variable. And since 

- (2) 

it is the volume-dependence of ip op t that renders P}-' 
negative, it must also be related to the "volume" terms 
discussed in Refs. @ ||, ||, || (see also Ref. JmJ). Our 
formulation of the problem hence illustrates that the way 
in which these terms drive a phase transition is closely 
related to the choice which variables one intends to keep 
constant when differentiating the thermodynamic poten- 
tial. Within the cell model it is obvious that PB theory 
only yields positive pressures and which pressure defi- 
nition is hence preferrable, but when relaxing the con- 
straints of a cell model or mean field PB theory the na- 
ture of effective interactions is much less obvious. See 
for instance Ref. j[l| for a recent critical evaluation of 
the theoretically as well as experimentally subtle issue of 
phase separation in suspensions of spherical colloids. 

In our next example we turn to a physical system 
which differs in the geometry of the colloids and which is 
(quite literally) of vital importance. Fig. [| shows pres- 
sure and compressibility of an aqueous solution of cylin- 
drical colloids, i.e., charged rods, which have a radius 
of ro = 1.2 nm and a line charge density of 1 charge per 
1.7 A (or 4.2 charges per Bjerrum length) dialyzed against 
an 1:1 electrolyte of concentration 100 mM. The colloid 
is much stronger charged than in the spherical situation 



above, in the sense that the surface charge density is a 
factor 50 larger. On the other hand, the salt concen- 
tration is also much larger (with the reservoir screening 
length being only 0.6% of the screening length in the 
above low-salt case) which keeps the potentials low. The 
above choice of values corresponds to DNA in a physio- 
logical salt environment. As Fig. |^ shows, nonlinear PB 
theory again gives positive pressure and compressibility 
for all volume fractions <j>, and the large-0 behavior is cap- 
tured correctly if Donnan linearization is employed, while 
the choice ib = fails there. However, using V'opt and the 
definition (|22|) results in a pressure that is negative for 
4> % 17% and a negative compressibility for <p < 11%. 
If this were true, all DNA in animal cells would tend to 
aggregate and phase separate! But again, this failure is 
not an inevitable artifact of linearization, since the pres- 
sure pfr^ is perfectly positive and gives rise to positive 
compressibilities. Incidentally, DNA can be condensed, 
but this requires multivalent ions |5lJ and is known to be 
a correlation effect which is missing in PB theory, see for 
instance || || |||. 

Poisson-Boltzmann theory and the linearization 
scheme employing -0 = predict an exponentially van- 
ishing pressure as <f> — ► 0. However, as we have remarked 
above the pressure due to the macroions has to become 
significant at some point (see Ref. J|5| for measurements 
on DNA in low salt solutions which appear to support 
this view). For linear polyelectrolytes this is even more 
relevant than for spherical colloids, since the conforma- 
tional degrees of freedom and the degree of entanglement 
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FIG. 4: Same plots as in Fig. |3|, but for cylindrical macroions having a charge parameter £ = 4.2 (i. e., 4.2 charges along the 
axis per Bjerrum length) and a radius ro = 1.2 nm (those values correspond roughly to DNA). The system is in contact with a 
dilute 1:1 electrolyte of physiological salt concentration 100 mM. The line styles are the same as in Fig. 0. 



becomes important | p6| . 

In our final example we study an asymmetric elec- 
trolyte. We go back to the system of spherical colloids 
studied above and replace the ions of one sign by diva- 
lent ones, i.e., we assume a 1:2 electrolyte of concen- 
tration 3.6 /iM, which in particular implies that the di- 
valent species occurs with a concentration of 1.8 p M. in 
the reservoir. According to the discussion in Sec. VI B 
the case in which the colloid has the same sign of charge 
as the divalent ions is particularly interesting, since then 
positivity of the pressure cannot be guaranteed even for 
the pressure definition ([21] ), For this case Fig. || shows 
again pressure and compressibility as a function of vol- 
ume fraction. PB theory is once more found to give pos- 
itive results, and the same remarks as above about the 
poor high-(/> behavior of the linearization point ip — 
compared to tp = tpopt a PPby also here. 

However, there is an (expected) difference concerning 
the behavior of P-f^ , which becomes (very slightly) neg- 
ative for <f) < 0.68% and gives rise to negative compress- 
ibilit y for < 0.51%, as is illustrated in the insets. In 
Sec. VI B we have seen that the pressure can become 
negative if \Q\ < 1.596nb^, which in this case implies 
4> ;$ 0.96% - and indeed, lowering <fi by further 30% brings 
about the negative pressure. We mention in passing that 
with a colloid having the opposite sign of charge as the 
divalent species this does not occur. 



VIII. CONCLUSION AND SUMMARY 

The main motivation of the present paper has been 
the following question: How can the membrane pressure 



due to a Donnan-equilibrium best be described within 
linearized PB theory? This first required a clarification 
of (i) how one computes the pressure and (ii) what one 
means by "best" . Often the pressure has been computed 
(without much justification) by using the predictions of 
the linearized PB equation in formulas of the nonlinear 
theory (or expansions thereof); moreover, this recourse to 
the full theory is not unique and hence prone to inconsis- 
tencies. These problems are avoided by starting from the 
appropriate thermodynamic potential functional of linear 
theory, which however still leaves two alternative defini- 
tions for the pressure that differ in their treatment of a 
possible volume dependence of the linearization point. 
Neither definition predicts the pressure to be propor- 
tional to the density of microions at the outer cell bound- 
ary, i. e., the well known result from nonlinear PB theory 
does not apply on the linearized level. All this can be dis- 
cussed without fixing the expansion point, but optimal- 
ity considerations nevertheless suggest a specific choice 
for it. Explicit formulas can be obtained, since the lin- 
earized PB equation can be solved analytically for general 
cell models. 

Let us thus sumarize the key findings of this work: 

(a) Linearized PB theory can be based on a density 
functional which is the quadratic expansion of the well 
known functional of nonlinear PB theory. The choice of 
the expansion point distinguishes different linearization 
schemes. The equilibrium value of this functional is the 
thermodynamic potential, in our case the grand poten- 
tial. 

(b) The pressure is given by the volume derivative of 
this thermodynamic potential. There is a choice as to 
whether or not one would like to keep the linearization 
point fixed, resulting in two different pressure definitions 
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FIG. 5: Same plots as in Fig. |3|, but for a 1:2 electrolyte. The colloid charge has the same sign as the divalent ions, the reservoir 
density of the monovalent ions is 3.6 fjM.. The inset in the left panel magnifies the regions [0; 0.008] x [—0.003; 0.003], while the 
inset on the right magnifies [10~ 4 ; 10~ 2 ] x [—0.01; 0.005], showing that in this example the first pressure definition ( |2l| ) also 
gives negative pressures and compressibilities, while the PB result does not. 



both based on the thermodynamic potential. 

(c) The boundary density rule @ from nonlinear PB 
theory is replaced by the quadratic expansion of the non- 
linear formula in the boundary potential. 

(d) The range of validity of linearization depends 
on the expansion point. We proved that linearization 
about the selfconsistently determined average potential is 
optimal — in the sense that it automatically describes the 
Donnan effect correctly in lowest order, leaving the lin- 
earized PB equation to incorporate higher order correc- 
tions. Furthermore, we showed that all other lineariza- 
tion schemes violate an important inequality from non- 
linear PB theory related directly to the Donnan effect. 

(e) The linearized PB equation can be solved exactly 
for symmetric cell models of arbitrary dimension, salt 
reservoir composition and linearization point. We used 
this solution to derive explicit formulas for the pressure 
and have discussed their analytical properties in detail. 

(f) A comparison with the results from nonlinear PB 
theory has been performed, showing that the validity of 
linearization depends strongly on the linearization point. 
While the traditional linearization scheme completely 
fails to describe the important limit of large volume 
fraction, the optimal (Donnan-) linearization becomes 
asymptotically correct there. 

(g) For pressure definition (^2|), i.e. the total volume 
derivative of the grand potential, negative pressures and 
compressibilities can occur. These findings have pre- 
viously been taken as indications of a gas-liquid phase 
separation in suspensions of spherical colloids. In the 
present work we showed that the same theoretical rea- 
soning would predict DNA to phase separate under phys- 
iological conditions at all relevant concentrations. We 



consider this as a further striking argument against the 
genuineness of the transition. 

(h) In nonlinear PB theory the pressure is always posi- 
tive. Whether or not this still holds on the linearized level 
depends crucially on the precise definition of the pressure 
itself. It is for this reason that we advocate defining the 
pressure via the partial derivative of the grand potential 
with respect to the volume (i.e., keeping in particular 
the linearization point fixed), since the pressure is then 
always positive for symmetric electrolytes — in agreement 
with nonlinear PB theory. 
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APPENDIX A: SOLUTION OF THE LINEAR PB 
EQUATION FOR A CELL MODEL IN d 
DIMENSIONS 

Consider a cell model in d dimensions with general- 
ized radial coordinate r. d = 1, 2, and 3 corresponds to 
planar, cylindrical and spherical macroions, respectively. 
After the transformation x = Rr and ip(r) = ip(Rr) = 
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"4>{x) the linearized PB equation (|15j) reads 

$"(x) + — i>'(x) = $(x) - Vi, (Al) 



where the prime indicates differentiation with respect to 
x. Henceforth we will not bother with the difference be- 
tween ip and ip an d om it the tilde. An obvious particular 
solution of Eqn. ( |Al[ ) is ip(x) = ipi, while the homoge- 
neous equation is solved by the ansatz if)(x) = w(x)/x v 
with v = d/2 — 1, provided w(x) solves the Bessel equa- 
tion x 2 w"(x) + x w'(x) — (x 2 + v 2 ) w(x) — 0. Two linearly 
independent solutions are the modified Bessel functions 
K„ and \ v Wh - The general solution of the differential 
equation (|Al|) can therefore be written as 



tp{x) = ipi + C K 



K v (x) 



(A2) 



where Ck and C\ are constants to be determined by the 
following boundary conditions. Gauss' law relates the 
radial derivative of the electrostatic potential at ro = 
xq /k to the surface charge density ea of the macroion. If 
we define the latter to be the total charge of the macroion 
divided by its surface area, we have ^'(^o) = —4tt£-b<j/k. 
Since the whole cell is neutral, we also have ip'(X) = at 
the outer cell radius R = X/R. The integration constants 
now follow from inserting these boundary conditions into 



the general solution. After a little algebra we then find 
the potential 



ip(x) = ipi + 



kD \ x 

l v+l {X)K v {x) +K v+l {X)l v {x)\, (A3) 

where the determinant D is defined as 

D = K v+1 {x Q )l v+1 {X)-K v+1 {X)l v+1 {x ). (A4) 

We have D > since X > x . 

At the outer boundary the potential (A3) can be fur- 
ther simplified. The term in angular brackets reduces to 
l/X, and using the volume fraction <fi the surface charge 
density can be rewritten in terms of the colloid charge as 
c/^o = (Q/V)(l — <j))/d(j). Inserting the definitions jl^ ) 
for R we arrive at 



ip(R) = Vi + 



Q/V 1 



(A5) 



Equation ( A5) gives the boundary potential in linearized 



PB theory for a general <i-dimensional cell model with 
arbitrary linearization point ip and electrolyte composi- 
tion. 
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